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We propose a new implementation of target mass corrections to nucleoli structure functions which, 
unlike existing treatments, has the correct kinematic threshold behavior at finite in the x ^ 1 
limit. We illustrate the differences between the new approach and existing prescriptions by consid- 
ering specific examples for the F2 and Fl structure functions, and discuss the broader implications 
of our results, which call into question the notion of universal parton distribution at finite . 
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I. INTRODUCTION 



\ Deep inelastic lepton scattering is one of the most developed tools with which to probe the quark and gluon structure 
of hadrons. The theoretical framework within which the experiments are analyzed is the operator product expansion 
' (OPE). This is well established in the Bjorken limit, which is defined by the four-momentum transfer squared and 
, energy transfer v being asymptotically large, with the ratio x — /2Mv (the Bjorken scaling variable) fixed, and M 
is the target mass. Within this framework, global analyses of deep inelastic scattering (DIS) and other experiments 
have allowed a vast array of data to be described in terms of a universal, process independent set of quark and gluon 
(or parton) distributions. 

On the other hand, there exists a large body of data at lower energies, at < 1 — 2 GeV^ (see e.g. Ref. ij), where 
. the use of the asymptotic, Bjorken limit formalism may be more questionable. In addition to perturbative QCD effects 
' generated by gluon radiation, at low effects arising from l/Q^ power corrections become increasingly important. 
I These are typically generated by multi-parton correlations, and within the twist expansion are associated with higher 
' twists (in the OPE the twist of an operator is defined as its dimension minus its spin). Since they characterize 
, the long-range nonperturbative interactions between quarks and gluons, the higher twists contain information on 
• confinement dynamics, and are as such of intrinsic interest to study. Several recent analyses of structure function 
"*7^ \ data at low have extracted matrix elements of the higher twist operators 2i] . 

'"^ . Before one can reliably extract information on the higher twist contributions, it is important to remove from the 
^ ' data corrections arising from purely kinematic effects associated with finite values of Q"^ jv^ = AM'^x^ /Q^ . These so- 
^ . called "target mass corrections" (TMCs) are formally related to twist-two operators, and hence contain no additional 
^ information on nonperturbative multi-parton correlations. Indeed, these have long been considered uninteresting and 

. ^ , believed to be well understood. In the literature there are well known prescriptions for how to remove the TMC 

^ ■ corrections 

\ The TMCs were first considered by Nachtmann 's^, who showed that one could arrange the OPE so as to ensure 
that at a given order in only operators of a given twist would appear. At finite the natural scaling variable, 

defined as the fraction of the nucleon's light-front momentum {p) carried by the parton [k), is then given by the 
Nachtmann scaling variable: 

^^^'^ = pOj^pz - 1 + ^i + 4M2a;2/Q2 ' 

In the Nachtmann approach, one generalizes the Cornwall-Norton (CN) moments of structure functions, derived in 
the Bjorken limit, to finite . A particular feature of these Nachtmann moments is that they are supposed to factor 
out the target mass dependence of the structure functions in a way such that its moments would equal the moments 
of the corresponding parton distributions. 

Later Georgi and Politzer (CP) jjj calculated the CN moments of the structure functions, taking into account the 
trace terms which appear in the matrix elements of the twist-two operators, but which are usually neglected in high 
data analyses. While the leading, twist-two piece of the structure function which enters at 0{1) is related to matrix 
elements of the quark bilinear ipj'^ip, for example, the TMCs arise from insertions of derivatives, ■07^I?^f'i • • ■ D~''^^ip, 
which does not alter the twist. This procedure generates a series in A'P/Q^ when calculating the CN moments. 
Inverting these moments using the inverse Mellin transform, one arrives at a structure function depending on both x 
and MVQ^ 0. 
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Problems with the GP implementation of TMCs were soon identified, however, by a number of authors |RF^IMI^ITol|. 
in particular the so-called "threshold problem" . This pertains to the fact that if the parton distribution function is a 
scaling function of ^, then since the maximum kinematic value of ^ at any finite is = ^(^^ = 1) < li the parton 
distribution is not defined in the unphysical region between the elastic limit ^ = and ^ = 1. De Rujula et al. [Tlj 
argued that the problems can be resolved by considering in addition higher twist operators. They note that there 
is a nonuniformity in the limits as rt ^ oo and — > cxd, and the appearance of higher twist effects proportional to 
nM^/Q^ for the n-th moment signals the breakdown of the entire approach at low (< 2 GeV). 

Tung and collaborators _8, 9] attempted to redress the threshold problem by invoking an ansatz which smoothly 
merges the perturbative QCD behavior of the moments at large with the correct threshold behavior in the n ^ oo 
limit. As they note, however, such a prescription is not unique, and in fact agrees with the standard OPE expansion 
only in the n — > cw limit. 

The proposed solution of De Rujula et al. to the threshold problem implies that higher twist effects play an 
important role at low Q^. Recent experiments at Jefferson Lab have shown, however, that the size of the higher twist 
effects is actually quite small for the proton F2 structure function, down to relatively low values (Q^ ^ 0.5 — 1 GeV^) 
p^. The question which we address here is whether a self-consistent formulation of TMCs can be made with only 
twist-two contributions, without appealing to higher twist effects. While not a proof, it seems plausible to us that, 
at least from a purely theoretical perspective, it should be possible to obtain an implementation of TMCs for a 
hypothetical case of negligible higher twist effects, which would demand a consistent resolution of the threshold 
problem for the twist-two part alone. Such a view could be motivated by observing that even though it is the same 
proton state that the twist-two and higher twist terms originate from, in principle the matrix elements of the local 
operators whose matrix elements characterize the twist expansion are in fact independent. 

While interesting in its own right, the question of how to implement TMCs is also of practical importance, given the 
high quality electron-nucleon structure function data at low and moderate which are being collected at Jefferson 
Lab hj. TMCs are also vital in analyzing neutrino scattering data , much of which are taken at relatively low 
energies, and must be understood if one is to extract reliable information on neutrino oscillations for instance. For 
spin-dependent scattering, TMCs have also been calculated for the gi and 32 structure functions , and recently for 
spin-1 targets such as the deuteron (141] . 

The pertinent question is whether the Nachtmann moment of the twist-two part of the structure function is 
independent, as supposed in the original formulation In Sec. II we review the standard derivation and results for 
TMCs within the operator product expansion. We outline the problems associated with the standard approach, and 
suggest an alternative formulation designed to avoid the unphysical threshold problem. Earlier work [lOj did indeed 
find that the Nachtmann moments do not account for all possible (leading twist) M^/Q^ effects. However, we suggest 
a prescription where the IVP/Q^ dependence of the Nachtmann moments of the structure functions does equal, to 
very high accuracy, the M^/Q^ dependence of the moments of the quark distributions, for all Q^. Numerical results 
are presented in Sec. Ill, where we compare the x dependence of the F2 and structure functions using the various 
TMC prescriptions, and examine the onset of scaling in terms of the Nachtmann moments of the structure functions. 
In Sec. IV we summarize our findings and discuss the broader implications of our results for the interpretation of 
parton distributions at finite Q^. 



II. OPERATOR PRODUCT EXPANSION 



We begin this section by firstly reviewing the pioneering work on target mass corrections as obtained by Georgi 
and Politzer We will consider the case of unpolarizcd scattering from a spin-1/2 nuclcon, which is described by 
two structure functions, Fi{x, Q^) and ^2(2;, Q^) (or alternatively F2 and the longitudinal structure function Fl). We 
shall focus on the F2 structure function, but later generalize the discussion to include also F^. 

The two standard moments of structure functions encountered in the literature are the Cornwall-Norton and Nacht- 
mann moments. The Cornwall-Norton moments of F2 are given by: 

M^{Q^) ^ r dx x"-^ F2{x,Q^) , (2) 
Jo 

and are appropriate for the region ^ Af^. The Nachtmann moments, on the other hand, take into account finite 
M'^/Q^ corrections to the Bjorken limit, and are given by: 

rl cn+l 

M^(g2) = dx^ 

Jo ^ 

with r = \/i + 4a;2M2/(52 ^j^g essential difference between the CN and Nachtmann moments comes from the trace 
terms appearing in the matrix elements of operators of definite spin, which are disregarded in the CN approach, but 



3 -h 3(n + l)r + n(?i + 2)H 
(n + 2)(n + 3) 



F2{x,Q^) , 



(3) 
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kept in the Nachtmann approach. The Nachtmann moments are constructed such that from the infinite operators of 
twist-2 and different spin contained in the trace terms, only the operators of spin n contribute for the n — 2 moment 
of the structure function. 

The Nachtmann and CN moments can be related by expanding the moments in powers of Expanding fj,2 to 

0{1/Q^), one has: 

,-.{Q') = M-iQ^) - ^^Q^Mr^Q^) + i^QT^^r'm - ^^QF^r^ + ••• (4) 

Note that there is a mixing between the lower and higher moments. To this order we can also express the CN moments 
in terms of the Nachtmann moments: 

M2{Q ) - ) + ^^^QjM2 {Q ) + 2(n + 3)(n + 4) Q^^' ^ 

n(n^-l)(n + 2)(n + 3) M6 „+g 
+ 6(n + 5)(n + 6) ^ "' ^ ' 

In the work of GP, the moment of the leading twist part of the F2 structure function, corrected for target mass 
effects, can be written to order as: 

^" ) = ^" + ^^Q^ + 2(. + 3)(U4)V ^""^^^^ 

n(n^-l)(n + 2)(n + 3) A/6 ^ 
+ 6(n + 5)(n + 6) Q6 + ■" ' 

where An is the n-th moment of a distribution function F[y): 

rya 

An = / ^^(y) . (7) 

Jo 

Here the function F{y) is related to the usual quark distribution q{y) by q{y) = yF{y), and the upper limit of 
integration j/o is the maximum value at which the quark distribution has physical support. Again in Eq. JBJ there is 
a mixture between lower and higher moments. Comparing Eqs. JSJ and 0, one can show that, at least to 0(1/ Q^), 
the Nachtmann moments are equivalent to the moments of the distribution F{y): 

^4=An. (8) 

This reflects the fact that the Nachtmann moments are constructed to protect the moments of the structure functions 
from target mass effects, thereby allowing them to be identified directly with the moments of the quark distributions. 

The F2 structure function appearing in Eqs. (0) and (jSjl must itself be corrected for target mass effects, and to this 
end we will follow the procedure in GP Q, albeit with one exception. While GP write the upper limit of integration 
in Eq. 10 as yo = 1; we will define the upper limit of the integrals as the maximum value allowed by kinematics, 
yo = y{x — 1). Following GP, we can then rewrite Eq. lO as Q: 

"^"^'^ - ^""dy G(y) , (9) 



(n + 2j)(n + 2j-l) 
with G{y) given by: 



fVo rVo rVo 

G{y) = / dy' H{y') = / dy' dy" F{y") . (10) 

Jy Jy Jyi 



This result follows from the fact that 

n+2j-2 



yo ... yn+2o-l fVo yn+2j-l QQ^y-^ 



and because G(0) = G{yo) = 0, one is left with the second term only. Integrating the RHS of Eq. again, we 
recover Eq. To obtain the x (or ^) dependence of the structure functions, we can invert the moment as in GP, 
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arriving at same results but with a modified upper limit in the integrals. Specifically, we recover for the F2 structure 
function 0]: 



6a 



HiO + 12a 



where a = M/Q. The CN moments of the target mass corrected F2 structure function are then given by 



dxx" 



2j 



j!(n-2)!(n + 2j)(n + 2j-l) 



To calculate the Nachtmann moments, we rewrite Eq. Q in terms of ^: 



.-2 (i + g^a^ 

1 - a2^2 



F2{LQ^) 



3(r-l) Or -If 



r2(n + 2) r2(n + 3) 



with i^2(C, Q^) given by Eq. (jTH), and 



(12) 



(13) 



(14) 



= C(x = 1) = 



1 + v/1 + 4M2/Q2 



(15) 



In the following section we will examine the extent to which the Nachtmann moments of F2 correspond to the moments 
An of the quark distribution function, for different functional forms of F(^), and quantify the effect of the kinematic 
thresholds. Central will be the interpretation of the function F(^) itself. 

Before proceeding, for completeness we also give the results for the longitudinal structure function, Fl, and its 
moments. In the Q2 — > 00 limit, F^ = 0, while at finite the TMCs render Fl nonzero. Of course, higher order 
perturbative QCD corrections which depend on as also give rise to a nonzero F^, as do higher twist effects. However, 
we artificially set both of these to zero in order to isolate the effects of TMCs on explicitly. Following a similar 
procedure as for F2 above, we can write the longitudinal structure function as: 



Fl{^,Q^) 



(1 + a2^2): 



(1 + ^2^2)3 



which clearly vanishes as a — > (or 00). The corresponding Nachtmann moments are then given by: 



Co 



4a2^2 („ + i)(i_a2^2)_2(n + 2) 



(l-a2^2) 



(n + 2)(n + 3) 



F2{^,Q' 



(16) 



(17) 



Having derived analytic expressions for the F2 and Fl structure functions and their moments, in the next section we 
present numerical results using several different prescriptions for the ^ dependence of the quark distribution function. 



III. TARGET MASS CORRECTIONS 



The main purpose of this work is to analyze phenomenologically what is the best procedure to incorporate TMCs 
in the analysis of structure functions. More specifically, we address the question of which procedure is most effective 
in rendering the moments of the leading twist structure functions equal to the moments of the quark distributions at 
finite (52 _ \Ye consider whether there is any sizable difference for the moments when the upper limit of the integrals 
in G{y) and H{y) is 1 or ^o- In particular, since the twist- two part of the deep inelastic cross section should be zero 
at a; = 1 (^ = ^o)j we study the impact of a vanishing parton distribution at ^o- 



A. Prescriptions 

To address these issues, in this section we present several prescriptions for the implementation of target mass 
corrections, and discuss their limitations and practical consequences. We consider three scenarios: 
(A) Integrate a quark distribution: 



qio = A/-r'/2(i-e)3 



(18) 
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from to 1 (specifically, in the integrals for An, H{(^) and G'(^)). Here the normalization Af ensures that the 
distribution integrates to unity. We denote this prescription the "standard TMC" (sTMC). 

(B) Integrate a modified distribution which vanishes for ^ > £,o, as implied by Eq. 0^: 

q{0 = Mr'^\i-0' e(c-^o) . (19) 

We denote this prescription the "modified TMC" (mTMC). 

(C) Use a "threshold dependent" (TD) quark distribution which vanishes in the physical limit: 

g™(0 = ^r'/'i^o-O'- (20) 



T ' 1 

n=2 



sTMC 
— mTMC 
TD 

_i 1 I 1 

12 3 4 

(GeV^) 

FIG. 1: Ratio of the n — 2 Nachtmann moment of the F2 structure function and the n = 2 moment of the quark distribution, 
as a function of Q^. The curves correspond to prescriptions A ["sTMC"] (dotted), B ["mTMC"] (dashed) and C ["TD"] (solid). 

Note that because of the upper limit in Eq. ||7J), An itself will be AP/Q^ dependent for prescriptions B and C. The 
results for the ratio 1^2/ An of the n = 2 moments are displayed in Fig. 1 for the three cases, with prescriptions A, B 
and C corresponding to the dotted, dashed and solid curves, respectively. Comparing the sTMC and mTMC results, 
one can see a reduced dependence when the integrals are restricted to ^ < ^o- However, a much more dramatic 
change occurs when the quark distribution is constrained to vanish at ^o- This renders the Nachtmann moment almost 
equal to the moment of the quark distribution for virtually all considered. Certainly for > 1 GeV^ there is no 
visible deviation of the ratio from unity. Even for very small Q^, ~ 0.3 GeV^, the ratio differs from unity by only 
^ 0.7% (of course the OPE itself may not be valid at such low values of Q^). 

Similarly, the ratios for the n — A and n — & moments are shown in Fig. [21 The deviation of the ratio from unity 
for the sTMC approach is between 10% - 20% for < 1 GeV^, while that for the modified TMC with prescription 
B is of the order of 5% over the same region. On the other hand, for the threshold dependent prescription C, the 
deviation from unity remains around 1% even at these low values. 

A consequence of prescription C is that the moments of the parton distribution are dependent. This seems to 
be an inevitable consequence if the Nachtmann moments of the structure function are to be equal to the moments of 
the parton distribution for all . Note that this dependence is not of higher twist or perturbative QCD origin, 
but arises solely from kinematics. Nevertheless, this avoids the more serious problems which arise within the sTMC 
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FIG. 2: Ratios of the n = 4 (upper graph) and n = 6 (lower graph) Nachtmann moment of the F2 structure function and the 
corresponding moments of the quark distribution, as a function of . The curves are as in Fig. 



approach (prescription A) , where the Nachtmann moments below Q'^ ^ 1 GeV^ start to deviate significantly from the 
moments of the quark distributions. In addition, in the sTMC formulation one is faced with the so-called "threshold 
problem". Namely, if the moments An of the quark distributions are independent, then one should have: 



Jo Jo 



) 



(21) 



for any two momentum scales Ql and Q2- Since F{^,Q^) must vanish in the kinematically forbidden region ^ > ^0, 
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the equality in Eq. (|21|l implies that the function must be zero for both ^ > Co(Qi) and ^ > ^0(^2)- 1^ Qi < Ql' i'^ 
which case Co(Qi) < ColQi); this implies that ^(CjQI) should vanish in the range £,o{Qi) < ^ < ^oiQ^)- However, 
there is no physical reason for it to do so here, and this therefore leads to an unphysical constraint. 

De Riijula et al. [Tl| address this problem by pointing out that higher twist contributions play an ever more 
important role at low Q^, and their neglect makes any leading twist analysis at large ^ incomplete. On the other 
hand, the philosophy inherent in prescription C is that the equality (|21|l should not be expected to hold, simply 
because the ^ dependence, and also the normalization, of the quark distributions is /Q^ dependent: 

?o(Q?) rioiQl) 

d^C F{^,Ql) + \ d^C F{^,Ql) . (22) 
Jo 

In the following section we will contrast the various prescriptions by studying their effects on the structure functions 
numerically. 

B. Numerical Results for Structure Functions 

The effects of TMCs on the F2 structure function are illustrated in Fig. O for — 1 and 5 GeV^ for the various 
scenarios. Here the scaling function xq{x) (dotted curve) is used as input in Eq. (|12|) to calculate the target mass 
corrected function F2. To translate the results from ^ to x, we fix and extract the corresponding x for each Q^) 
pair. For the sTMC method (prescription A, dashed), in which the upper limits of the integrals are set to unity, the 
corrected structure function becomes smaller at intermediate x values, but larger as a; — > 1. The corrected structure 
function for the mTMC approach (prescription B, double-dot-dashed), where the ^ integration is constrained by 
C < display a similar behavior as a function of x. In both cases F2 is clearly nonzero in the 2: — > 1 limit. The 
effects are sizable at = 1 GeV^, but considerably smaller a,t — 5 GeV^, where the differences between the 
scaling and target mass corrected functions are more strongly suppressed. At lower values, 0.5 GeV^ (not 

shown), the differences between the sTMC and mTMC prescriptions are even more pronounced, so that here it is 
important to take into account the correct kinematics, especially at large x. Note that the sTMC and mTMC curves 
in Fig. 121 were normalized such that the quark number at finite is equal to the quark number at 00. 

The effect of the threshold dependent prescription C is dramatically different from the other prescriptions. Specif- 
ically, the TD input distribution q^^ (dot-dashed curve in Fig. |2Jl produces a target mass corrected leading twist 
structure function which is exactly zero a,t x = 1 (C = Co) for all Q^, as required physically. Another important 
difference between the sTMC and mTMC approaches, and the TD prescription, is at intermediate x. Here the 
latter produces a corrected structure function which is smaller than the input, in contrast with the sTMC/mTMC 
methods, where the corrected F2 is larger than the input scaling function. Such differences may be very relevant in 
phenomenological determinations of the x dependence of parton distributions at low Q^. 

How does the TD prescription C help in the practical extraction of leading twist parton distributions? If we define: 

A„^Ccr("+^/^)A.(Q'), (23) 

then since at high the moments of the quark distributions become independent, so too are the Nachtmann 
moments of the structure functions, as the multiplicative factor is C independent. Other prescriptions would render 
a different Co: but the important fact is that once (Q^) = A„(Q^) for any Q^, this equality may be set to its value 
at — * 00: 

/x^ (finite g^) ^ /ig(Q2 ^ 00) ^ Ml 

^„ (finite g2) A„(Q2^oo) A„ ' ^ ' 

This makes the approach in prescription C much more useful for extracting quark distributions from structure function 
data at low Q^. 

A comparison of the different TMC prescriptions for the longitudinal structure function moments is shown in Fig. 01 
for the n = 2 and n = 4 Nachtmann moments. Here we plot the absolute value of the moments rather than a ratio, 
since the moments in the scaling limit are identically zero. For the sTMC prescription, the corrected moments are 
strongly dependent for ^ 1 GeV^ and rise rapidly as decreases. The dependence of the moments 
for the mTMC approach, with the upper integration limit being Co rather than unity, is somewhat weaker, but still 
quite strong at low Q^. On the other hand, the moments for the TD prescription C display significantly smaller 
dependence at low Q^. 

Finally, for completion, we show in Fig. [S] the structure function at = 1 GeV^ with TMCs applied according 
to the three prescriptions above. Note that in the scaling limit, the structure function is identically zero. For 
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FIG. 3: The x dependence of the F2 structure function at = 1 GeV^ (upper) and 5 GeV^ (lower). The effects of TMCs 
on the (input) scaling distribution (dotted curve) are illustrated for the sTMC (dashed) and mTMC (double-dot-dashed) 
prescriptions, and compared with the effects on the (input) TD-distribution £,q^^{^) (dot-dashed) using the TD approach 
(prescription C, solid). 

the sTMC and mTMC prescriptions, the corrected structure function is significantly larger in magnitude than for the 

TD prescription at intermediate and large x. For the sTMC case in particular, it is also seen to approach a nonzero 
value in the x ^ 1 limit. This result suggests that the evaluation of the twist-two part of the longitudinal structure 
function at low may also need to be reassessed in phenomenological analyses, especially at intermediate and large 

X. 
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FIG. 4: Nachtmann moments of the longitudinal structure function for n = 2 (upper) and n = 4 (lower) as a function of 



IV. CONCLUSION 

In this work wc have revisited the long-standing problem of target mass corrections to nucleon striicture functions. 
The standard procedure for implementing target mass effects suffers from the well known threshold problem, in 
which the corrected, leading twist structure function does not vanish at x = 1. We have proposed a solution to this 
problem by introducing a finite-Q^, "threshold dependent" parton distribution function that explicitly depends on the 
kinematical threshold ; which is smooth in the entire physical region, and approaches the ordinary, Q^-independent 
parton distribution in the limit — > oo. Our prescription avoids any discontinuities in the parton distributions and 
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FIG. 5: Longitudinal structure function Fl at = 1 GeV^ for the sTMC (dashed), mTMC (double-dot-dashed) and TD- 
distribution (solid) prescriptions. Note that the scaling longitudinal distribution is zero. 

structure functions at finite Q^, and produces vanishing structure functions as a; — > 1. This is true for both the F2 
and Fl structure fmictions. 

The Nachtniann moments /ij of the F2 structure function, calculated with the threshold dependent distributions 
q"^^, agree with the moments An of q^^ to within 1% for the n = 2,4 and 6 moments for as low as 1 GeV^ and 
even lower. In contrast, the deviation for the standard or modified TMC procedure (sTMC or mTMC prescriptions) 
is more than an order of magnitude larger at the same values, and grows rapidly with increasing n. Furthermore, 
for > one can show analytically that, at least to 0{1/Q^), the moments ^^'^ are identical. Similarly, 
for the longitudinal structure function Fl, the Nachtmann moments /xj with the threshold dependent distribution 
are considerably smaller {i.e. closer to the asymptotic value of zero) than the moments in the sTMC or mTMC 
prescriptions. 

A consequence of our formulation is that the moments of the threshold dependent distributions will in general be 
AP/Q^ dependent. This dependence is not associated with either perturbative QCD effects or higher twists, but 
comes entirely from the leading twist, target mass effects. Our analysis suggests that it may be necessary to reassess 
the interpretation of a parton distribution in the presence of the finite M^/Q^, or ^, corrections, as well as the 
imp lementation of the q^^ distributions in the evolution equations. We will address these problems in future work 
UJ. At the same time, our numerical results give impetus to investigating the impact of TMCs on phcnomcnological 
fits to structure functions at low 'l6'| and the extraction of twist-two parton distributions. 
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